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ABSTRACT 
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Reynolds  stress  calculations  of  homogeneous  turbulent  shear  flow  are  conducted  with  a 
second-order  closure  model  modified  to  account  for  non-equilibrium  vortex  stretching  in  the 
dissipation  rate  transport  equation  as  recently  proposed  by  Bernard  and  Speziale  [./.  Fluids 
Fiif/rif/.  114,  29  ( 1992)].  As  with  the  earlier  reported  l\  —e  model  calculations  incorporating 
this  vortex  stretching  effect,,  a  product ion-equals-dissipation  equilibrium  is  obtained  with 
hounded  turbulent  kinetic  energy  and  dissipation.  However,  this  equilibrium  is  now  not 
achieved  until  the  dimensionless  time  Si  >  60  an  ('lapsed  time  that  is  at  least  twice  as  large 
as  any  of  those  considered  in  previous  numerical  and  physical  experiments  on  homogeneous 
shear  flow.  Direct  quantitative  comparisons  between  the  model  predictions  and  the  results 
of  experiments  arc  quite  favorable.  In  particular,  it  is  shown  that  the  inclusion  of  this  non¬ 
equilibrium  vortex  stretching  effect  has  the  capability  of  explaining  the  significant  range  of 
production  to  dissipation  ratios  observed  in  experiments. 


•Tin  s  research  was  supported  l>y  the  National  Aeronautics  and  Spare  Administration  under  NASA  Con- 
t  ract  No.  NAS  I  - 1  StJOT)  and  NASI- 1!)  ISO  while  the  first  author  was  in  residence  at  I  he  Inst  it  ute  for  ( ’omputer 
Applications  in  Science  and  lOngiiieering  (ICASK),  NASA  Langley  Research  Center,  Hampton,  VA  'idtitio . 
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INTRODUCTION 


In  transport  models  for  the  turbulent  dissipation  rate  it  is  traditionally  assumed  that,  an 
equilibrium  exists  where  the  production  of  dissipation  by  vortex  stretching  is  exactly  counter¬ 
balanced  by  the  leading  order  part  of  the  destruction  of  dissipation  term  (see  Tennekes  and 
Lumley  1972  and  Lesieur  1990).  Recently,  Bernard  and  Speziale  (1992)  proposed  a  new 
transport  model  for  the  dissipation  rate  with  residual  vortex  stretching  to  account  for  small 
departures  from  equilibrium.  This  new  model  -  which  is  of  the  same  general  form  as  that 
obtained  from  the  theory  of  self-preservation  (see  Speziale  and  Bernard  1992  and  George 
1992)  -  was  shown  to  yield  a  more  complete  description  of  isotropic  decay.  More  importantly, 
however,  it  was  demonstrated  by  Bernard  and  Speziale  (1992)  that  the  inclusion  of  this 
vortex  stretching  effect  yields  a  production-equals-dissipation  equilibrium  in  homogeneous 
shear  flow,  with  bounded  turbulent  kinetic  energy  and  dissipation.  Prior  to  saturation  - 
which  occurs  at  an  elapsed  time  larger  than  any  of  those  considered  in  previously  conducted 
physical  and  numerical  experiments  -  the  turbulent  kinetic  energy  and  dissipation  rate  grow 
exponentially  with  time.  Good  qualitative  agreement  between  the  model  predictions  with 
vortex  stretching,  and  the  results  of  physical  and  numerical  experiments  on  homogeneous 
shear  flow,  was  demonstrated  by  Bernard  and  Speziale  (1992)  and  Bernard,  Thangam  and 
Speziale  ( 1992).  However,  they  were  not  able  to  make  more  detailed  quantitative  comparisons 
with  experiments  since,  for  simplicity,  they  based  their  calculations  on  the  K  —  s  model.  This 
establishes  the  motivation  of  the  present  paper:  to  present  full  Reynolds  stress  calculations  of 
homogeneous  shear  flow  with  a  state  of  the  art  second-order  closure  model  that  incorporates 
this  vortex  stretching  effect. 

In  this  paper,  more  detailed  calculations  of  homogeneous  shear  flow  are  presented  based 
on  the  second-order  closure  model  of  Speziale,  Sarkar  and  Gatski  (1991)  modified  to  account 
for  non-equilibrium  vortex  stretching  in  the  dissipation  rate  transport  equation.  For  elapsed 
times  Si  <  30,  vortex  stretching  has  little  qualitative  effect  on  the  solut  ion  which  exhibits  an 
exponential  time  growth  of  turbulent  kinetic  energy  and  dissipation  that  is  in  good  agreement 
with  the  results  of  physical  and  numerical  experiments  (see  Tavoularis  and  Karnik  1989  and 
Rogers,  Moin  and  Reynolds  1986).  However,  for  St  >  60  (which  is  far  beyond  the  largest 
elapsed  time  that  has  been  considered  in  experiments),  the  flow  saturates  to  a  production- 
equals-dissipation  equilibrium.  Interesting  features  of  this  vortex  st  retching  solut  ion  which 
include  its  apparent  ability  to  explain  the  wide  range  of  production  to  dissipation  ratios 
observed  in  experiments  -  are  discussed  along  with  the  implications  for  turbulence  modeling. 


FORMULATION  OF  THE  PROBLEM 


We  consider  homogeneous  turbulent  shear  flow  where  the  mean  velocity  gradient  tensor 
takes  the  form 

P-  =  ■ifi.fn  in 

().Vj 

given  that  .S'  is  a  constant  shear  rate  and  htJ  is  the  Kronecker  delta.  In  the  usual  formulations 
of  homogeneous  shear  flow,  an  initially  decaying  isotropic  turbulence  is  subjected  to  the 
uniform  shear  rate  S  at  time  /  =  0.  The  Reynolds  stress  tensor  rtJ  =  u^u'j  is  a  solution  of 
the  transport  equation  (c.f.  Hinze  1975) 

()Tij  ()u,  2  c 

T“  =  + "« “  r(,i  (2) 


in  incompressible  homogeneous  turbulence,  where  a  superposed  dot  represents  a  time  deriva¬ 


tive  and 


dll';  (hl'i 


\d*i  d.r,r  -dxjdxj 

denote,  respectively,  the  pressure-strain  correlation  and  the  turbulent  dissipation  rate.  In 
(•'?).  //  and  u't  represent  the  fluctuating  pressure  and  fluctuating  velocity,  respectively,  and  // 
denotes  the  kinematic  viscosity  of  the  fluid:  as  with  most  model  studies  of  high-Reynolds- 
uumber  turbulence,  Kolmogorov's  assumption  of  local  isotropy  is  invoked  for  the  dissipation 
rate. 

In  order  to  achieve  closure,  models  must  be  provided  for  the  pressure-strain  correlation 
FI,,  and  the  turbulent  dissipation  rate  f.  The  recent  pressure-strain  model  of  Speziale,  Sarkar 
and  (lat.ski  (1991)  is  utilized  which  takes  the  form: 
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given  that  btJ  is  the  anisotropy  tensor,  A  is  the  turbulent  kinetic  energy,  and  V  is  the 
turbulence  production.  The  constants  of  the  model  are  given  as  follows:  (\  —  3.4,  fj*  =  1.80. 
("2  =  4.2,  (  3  =  4/5,  f’j  =  1.30.  ( ',t  =  1.25  and  (\  =  0.40.  Equation  (2)  is  solved  with  the 
modeled  dissipation  rate  transport  equation 

f  =  Cel ~v  +  (8) 

A  A  A 

where  Rt  =  K 2 / vs  is  the  turbulence  Reynolds  number  and  (\\ .  (\;  and  ( ',.-3  are  constants 
(see  Bernard  and  Speziale  1992).  The  term  containing  Oe3  occurs  when  there  is  a  departure 
from  equilibrium  so  that  there  is  an  imbalance  between  the  production  of  dissipation  by 
vortex  stretching  and  the  leading  order  part  of  the  destruction  of  dissipation  term  which 
each  scale  as  Rl/\  Equation  (8)  is  of  the  same  general  mathematical  form  as  that  obtained 
from  the  theory  of  self-preservation  (see  Speziale  and  Bernard  1992  and  Cleorge  1992).  The 
standard  modeled  dissipation  rate  transport  equation  is  recovered  in  the  limit  as  Ce: *  goes  to 
zero.  The  same  values  of  and  ( V2  as  proposed  in  the  Speziale,  Sarkar  and  (latski  (SS(!) 
model  are  chosen: 

Cel  =  1.44.  =  1 .83.  (9) 

A  variety  of  values  of  (\,\  in  the  range  of  0.001  to  0.01  will  be  considered  (these  represent t 
small  imbalances  in  vortex  stretching  of  the  order  of  0.1  %  -  1.0%;  of  course,  the  magnitude 
of  this  imbalance  can  depend  on  the  initial  conditions). 

The  governing  nonlinear  differential  equations  for  homogeneous  shear  flow  are  obtained 
by  substituting  (1 )  into  (2)  and  (8)  while  making  use  of  the  pressure-strain  model  (4).  Five 
coupled  nonlinear  differential  equations  for  ru ,  192,  r-n.  r;? 3,  and  £  comprise  t  his  system.  Those 
coupled  equations  can  be  easily  converted  to  an  equivalent  set  in  terms  of  h\  \ ,  b\2-  ^22 •  ^33 
and  Sl\/s  which  constitute  the  non-dimensional  structural  parameters  of  the  problem  that 
achieve  equilibrium  values  that  are  independent  of  the  initial  conditions  (see  Speziale  and 
Mac  (iiolla  Mhuiris  1989a,b).  The  governing  system  of  differential  equations  are  solved 
subject  to  the  initial  conditions: 

S  A  S  A  n 

K  =  0.  Rt  =  Rt0,  - =  — -  (,0) 

t  CO 

at  time  /  =  0.  Only  the  solutions  with  net  vortex  stretching  ((  Vi  >  0)  require  the  specifica¬ 
tion  of  the  initial  turbulence  Reynolds  number. 


DISCUSSION  OF  THE  RESULTS 


We  now  present,  results  obtained  from  a  Runge-Kutta  numerical  integration  scheme.  In 
Figure  1(a).  the  time  evolution  of  the  dimensionless  turbulent  kinetic  energy  (A*  =  A’/A’o) 


is  shown  as  a  function  of  the*  dimensionless  time  (/*  =  St)  for  SI\a/c0  =  3.38  and  Rlo  =  750 
the  initial  conditions  of  the  large-eddy  simulation  of  Bardina,  Ferziger  and  Reynolds  ( 1983). 
Here,  we  set  (  S}  =  0.001  so  that  there  is  only  a  minute  imbalance  in  vortex  stretching  of  the 
order  of  0.1fX.  As  with  the  earlier  reported  calculations  with  the  SS(1  model  where  ( ',3  =  0, 
there  is  excellent  agreement  between  the  model  predictions  and  the  large-eddy  simulation 
results  as  illustrated  in  Figure  1(a).  However,  while  there  is  virtually  no  difference  between 
the  short-time  predictions  of  the  SSC!  model  with  Ce 3  =  0  and  Ce 3  =  0.001,  there  is  a  distinct 
difference  in  the  long-time  solutions.  When  Ce 3  =  0.  the  kinetic  energy  and  dissipation  rate 
grow  unbounded  with  time;  for  /*  1, 

A'*  ~  rx,\  em  ~  eA*‘  (11) 

so  that  A*  and  c*  — >  oc  as  /*  — >  00.  For  any  finite  Ce 3  -  no  matter  how  small  it  may  be  - 
the  turbulent  kinetic  energy  and  dissipation  rate  eventually  saturate  to  bounded  equilibrium 
values  after  an  early  time  exponential  growth  as  shown  in  Figure  1(b). 

In  Figure  2.  the  time  evolution  of  the  turbulent  kinetic  energy  is  shown  for  the  same 
initial  conditions  (Sh'o/zo  =  3.38,  Rto  =  750)  but  with  a  variety  of  different  values  for 
("fi.  It  is  clear  from  these  results  that  the  value  of  Cc, 3  determines  the  saturation  level  of 
the  turbulent  kinetic  energy.  For  sufficiently  small  Ce 3,  the  turbulent  kinetic  energy  grows 
exponentially  for  St  <  30  the  largest  elapsed  time  considered  in  all  previously  conducted 
physical  and  numerical  experiments.  Here,  saturation  does  not  occur  until  St  >  60,  which 
is  a  value  that  is  more  than  double  that  of  any  considered  in  previous  experiments.  If  the 
vortex  stretching  solution  were  examined  in  isolation  for  elapsed  times  St  <  30.  it  could 
be  erroneously  concluded  that  there  is  an  unbounded  exponential  time  growth  of  turbulent 
kinet  ic  energy.  This  points  to  the  danger  of  drawing  conclusions  about  long-time  asymptotic 
states  from  data  corresponding  to  only  limited  elapsed  times. 

In  Figure  3(a).  the  time  evolution  of  the  ratio  of  production  to  dissipation  V  /  £  predicted 
by  the  SS(  1  model  both  with  and  without  vortex  stretching  is  shown.  The  non-zero  values  of 
the  constant  ( ',3  are  in  the  range  of  0.001-0.009  (an  imbalance  in  vortex  stretching  of  the  order 
of  0.1f7M(7f).  It  is  clear  from  this  figure  that  just  a  minute  imbalance  in  vortex  stretching 
can  cause  a  significant  spread  in  V /e  at  St  =  30.  This  is  reminiscent  of  the  experiments  of 
Tavoularis  and  Karnik  (1989)  where,  for  comparably  large  shear  rates  S l\  0/-0  >  1.  values 
of  V / £  in  the  range  of  1.1  -  1.8  have  been  observed  at  the  latest  station  measured  which 
corresponds  to  St  <  30.  Furthermore,  V/s  peaks  at  St  ~  20  before  asymptoting  to  a  value  of 
P/-:  =  1  which  is  achieved  for  values  of  St  >  60  as  shown  in  Figure  3(b).  This  phenomenon  of 
a  local  maximum  being  reached,  followed  a  gradual  dropoff  in  the  time  interval  0  <  St  <  30. 
is  reminiscent  of  results  observed  in  physical  and  numerical  experiments  (see  the  direct 
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simulations  of  Rogers  et  al.  1986  shown  in  Figure  4).  Without  this  vortex  stretching  effect, 
second-order  closures  predict  that  V/e  goes  to  an  equilibrium  val  ueof  (a*-  i)/(reI  -  1) 
monotonically ,  getting  close  by  the  time  St  —  20. 

As  in  the  absence  of  this  vortex  stretching  effect,  second-order  closure  models  predict 
that  btJ  and  S  K  /  s  achieve  equilibrium  values  that  are  independent  of  the  initial  conditions 
as  well  as  the  constants  C7e,,  Cc 2  and  CeJ.  In  Table  1,  the  equilibrium  values  predicted 
by  the  standard  SSCI  model  (Ce 3  =  0)  and  the  SSCI  model  with  vortex  stretching  ( > 
0)  are  compared  with  the  physical  and  numerical  experiments  of  Tavoularis  and  Karnik 
(1989)  and  Rogers  et  al.  (1986).  The  inclusion  of  vortex  stretching  only  causes  modest 
deviations  from  the  experimental  values  for  the  anisotropies;  only  the  shear  parameter  is 
affected  substantially.  However,  it  must  be  remembered  that  the  physical  and  numerical 
experiments  have  only  been  conducted  for  St  <  30.  For  the  case  where  Sh’o/co  =  3.38, 
R,0  =  750  and  Ce%  =  0.001  shown  in  Figure  1,  the  SSCI  model  with  vortex  stretching 
predicts  that 

/>,,  =  0.215,  fc12  -  -0.163,  bn  =  -0.141,  633  =  -0.074.  SK/e  =  5.40  (12) 

at  St  —  .30  results  that  are  in  close  proximity  to  the  experimental  data  which  is  given  011 
Table  1. 

CONCLUDING  REMARKS 

The  Reynolds  stress  calculations  presented  in  this  study  clearly  indicate  that  the  in¬ 
troduction  of  the  non-equilibrium  vortex  stretching  effect  of  Bernard  and  Speziale  (1992) 
in  a  recently  proposed  second-order  closure  model  does  not  compromise  the  accuracy  of  its 
quantitative  predictions  for  homogeneous  shear  flow.  In  fact,  when  the  SSCI  model  was  mod¬ 
ified  to  incorporate  this  vortex  stretching  effect,  better  agreement  with  experimental  data 
was  achieved  on  two  fronts:  (a)  consistent  with  experiments,  a  universal  equilibrium  was 
not  reached  by  St  =  30,  and  (b)  a  variety  of  production  to  dissipation  ratios  in  the  range 
1  <  V/e  <  1.8,  depending  on  the  initial  conditions,  were  obtained  at  St  =  30  consistent 
with  experiments  {V/e  asymptotes  to  1  for  elapsed  times  St  >  60).  In  our  opinion,  the 
level  of  agreement  between  the  model  predictions  and  experiments  is  of  such  quality  that  a 
product  ion-equals-dissipat  ion  equilibrium,  with  bounded  energy  states,  cannot  be  ruled  out 
as  a  viable  alternative  to  the  classically  accepted  hypothesis  of  unbounded  energy  growth. 

As  discussed  by  Bernard  and  Speziale  ( 1992),  the  question  of  whether  or  not  homogeneous 
shear  flow  actually  saturates  to  a  production-equals-dissipation  equilibrium  remains  an  open 
question  that  will  only  be  firmly  resolved  by  a  rigorous  mathematical  analysis.  However,  even 
if  it  were  to  be  ultimately  proved  that  homogeneous  shear  flow  has  unbounded  energy  growth. 
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this  would  l)o  duo  to  the  idealized  nature  of  t  ho  problem;  real  physical  systems  saturate 
before  singularities  occur.  The  introduction  of  this  vortex  stretching  effect  in  Reynolds 
stress  closures  could  lead  to  better  behaved  models  by  preventing  such  singularities,  without 
compromising  their  ability  to  collapse  experimental  data  in  benchmark  turbulent  flows  like 
homogeneous  shear  flow. 
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Equilibrium 
Val  ues 


SSC  Model 
(Cc3  =  0) 


SSC  Model 

(C,3  >  0) 


Experimental 

Data 


bn 

0.219 

0.201 

0.21 

bn 

-0.146 

-0.127 

-0.14 

h  3 

-0.073 

-0.074 

-0.07 

bn 

-0.164 

-0.160 

-0.16 
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S  I\  ft 
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5.0 

Table  1.  Comparison  of  the  predictions  of  the  SSC  model  for  the  equilibrium  values  in 
homogeneous  shear  flow  with  physical  and  numerical  experiments  (Tavoularis  and  Karnik 
1989  and  Rogers  et  al.  1986). 
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figure  1.  Tune  evolution  of  the  turbulent  kinetic  energy  in  homogeneous  shear  (low:  ,S  h0/£0 
=  3.38  and  /?<„  =  750.  (a)  Comparison  of  the  vortox-strctching  modified  SSC  model  ( ( ’,3  = 
0.001)  with  the  large-eddy  simulation  of  Bardina  et  al.  (1983),  and  (b)  the  long-time  SSC 
model  predictions  (Ce 3  =  0  and  Ce 3  =  0.001). 
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I1  igure  1.  I  ime  evolution  of  the  turbulent  kinetic  energy  in  homogeneous  shear  flow:  Sho/^o 
=  3-38  and  Rto  =  750.  (a)  Comparison  of  the  vortex-stretching  modified  SSCJ  model  (Ce^  = 
0.001)  with  the  large-eddy  simulation  of  Bardina  et  al.  (1983),  and  (b)  the  long-time  SSCI 
model  predictions  (Ce3  =  0  and  (7e3  —  0.001). 
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Figure  3.  J  line  evolution  of  tlie  ratio  of  production  to  dissipation  in  homogeneous  shear  flow 
for  S  A0/e0  =  3.38  and  Rto  =  750:  (a)  Predictions  of  the  SSG  model  for  a  variety  of  values 
of  for  St  <  30,  and  (b)  Predictions  of  the  SSG  model  for  St  <  100. 
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I'  igl,rr‘  I-  I  hue  evolution  of  Uio  ratio  of  production  to  dissipation  in  homogeneous  shear  (low 
for  Sl\0/£r>  =  .‘5.38  and  Rto  =  750:  (a)  Predictions  of  the  SSG  model  for  a  variety  of  values 
of  75.1  for  St  <  .30,  and  (l>)  Predictions  of  the  SSG  model  for  St  <  100. 
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